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Abstract 

The paper deals with optical excitations arising in a one-dimensional chain of identical spheres 
due optical coupling of whispering gallery modes (WGM). The band structure of these excitations 
depends significantly on the inter-mixing between WGMs characterized by different values of an- 
gular quantum number, I. We develop a general theory of the photonic band structure of these 
excitations taking these effects into account and applied it to several cases of recent experimental 
interest. In the case of bands originating from WQMs with the angular quantum number of the 
same parity, the calculated dispersion laws are in good qualitative agreement with recent experi- 
ment results. Bands resulting from hybridization of excitations resulting from whispering gallery 
modes with different parity of I exhibits anomalous dispersion properties characterized by a gap in 
the allowed values of wave numbers and divergence of group velocity. 

PACS numbers: 42.60.Da, 42.82.Et, 42.70.Qs 
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I. INTRODUCTION 



Recent proposal of coupled resonator optical waveguides and optical filters 1 ' 2 stimulated 
interest in systems of optically coupled micro-spheres. It has been known for a long time that 
electromagnetic modes of an individual sphere (Mie resonances) with large enough values 
of their angular momentum can have very long radiative life-times^. These are so called 
whispering gallery modes (WHM) characterized by the concentration of the electromagnetic 
field along the surface of a sphere with an evanescent tail escaping outside. Such field 
configuration makes it possible to optically couple two spheres positioned in the proximity 
of one another. Initial work on the optical coupling of the micro-spheres was concentrated 
on the case of just two spheres (the arrangement is known as photonic atoms or molecules), 
where the splitting of the modes and the formation of the coupled states was observed^*^. 
Recently, however, an interest has shifted toward linear chain of many spheres, which are 
envisioned as building blocks of various photonic devices, such as waveguides 7 or delay 
lines 8 -. The obvious result of optical coupling in this system is formation of collective optical 
excitations propagating along the chain, which we will call supermodes in order to distinguish 
them from the modes of individual spheres. The supermodes form photonic bands, which 
were observed by several research groups almost at the same tim o 7 i 8 i 9 . 

One of the popular tools used to analyze experimental results 8 is a simple phenomeno- 
logical dispersion law of a tight-binding type 

u = uq + k cos(qd) (1) 

where u and Uq are respectively the frequency of the supermode, characterized by a wave 
number q, and that of a single mode WGM resonance; phenomenological coefficient k char- 
acterizes the strength of the optical coupling, and d is the period of the structure. This 
approach, which was originally suggested in Refill to describe modes of coupled cavities, 
derive every photonic band from a respective single sphere WGM resonance, characterized 
by its angular momentum, I. Its applicability is based on the assumption that the modes of 
coupled spheres with different / do not mix, which is certainly not true in the exact sense, 
but can be approximately valid or not valid at all depending upon the type of modes under 
consideration. 

The admixture of WGM with different angular momentums arises already for two in- 
teracting spheres, and results from the violation of spherical symmetry in such systems. 
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It has been realized early on that the inter-mode coupling is primary responsible for ra- 
diative decay of the coupled modes- , and may also affect positions of the resonances in a 
bi-spherical structure. We will show in this paper that these effects are even more impor- 
tant, and in certain situations, crucial, for the supermodes of the multi-spherical chains. In 
these systems, instead of individual resonances, one has to consider mixing between bands 
of collective propagating excitations. Using analogy with solid state physics one can call this 
phenomenon band-mixing or band hybridization. In this paper, we present a theory of the 
inter-band coupling in linear chain of spheres based on the tight-binding approach to the 
optical coupling. This approach was originally formulated in Refllll. and carefully analyzed 
for the case of bi-spheres in Refill. 

Our main goal in this paper is to extend the analysis of Refill to a linear chain of 
spheres, and apply it to the problem of photonic band structure and dispersion laws of the 
quasi- stationary (long-living) collective optical excitations of the chain. Using a perturbative 
approach, we derive analytically dispersion equations characterizing these excitations and 
analyze them in several particular cases of recent experimental interest. We show, in partic- 
ular, that in spite of a formal analogy between optical and electronic tight-binding models, 
the photonic band mixing may lead to very unusual effects, which do not have counterparts 
in solid state systems. 



II. TIGHT-BINDING MODEL FOR A LINEAR CHAIN OF DIELECTRIC 
SPHERES 

The system, which will be considered in the paper consists of a number of identical spheres 
with radius R and refractive index n, whose centers are all aligned along the same line (see 
Fig^) at a distance d > 2R from each other. We will be interested in collective excitations 
of this chain, which in the spirit of tight-binding approach is described as a combination of 
modes of individual spheres. It is pertinent, therefore, to start with a brief review of the 
properties of a single sphere. Electromagnetic field in this case can be described by vector 
spherical harmonics (VSH), whose angular dependence is specified by angular number I 
(angular momentum) and azimuthal number m (z-component of the angular momentum). 
The radial part of the VSH, fi(kr) is given either by spherical Hankel function of the first 
kind, h\ (kr), for points outside of the sphere or by spherical Bessel function, ji(nkr), for 
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FIG. 1: The chain of spheres of radius R with Z-&xis along the axis of the chain 



points inside the sphere. Parameter k here represents vacuum wave number, k — u/c, where 
to and c are frequency and the speed of light in vacuum. We will use notation N m ^(r) and 
M m) /(r) for VSH describing electromagnetic waves of TE and TM polarizations respectively. 
Explicit expressions for functions N m /(r) and M m ;(r) can be found in many sources (see, 
for LJ! RefQ thns we do not pin, it here.' An i 8ol ated di.ect.c JL - an open 
system, therefore, it is more prudent to describe its properties in terms of resonances rather 
than normal modes. To this end one presents an incident, Ei nc , and scattered by the sphere, 
E s , waves in terms of VSH: 



E inc = [Cl,m^m,l( r ) + Vl,nMm,l( r )} ( 2 ) 

l,m 

E s = J2 h,™N m>i (r) + 6/, m M m ,i(r)] . (3) 

l,m 

Using Maxwell boundary conditions one establishes relations between coefficients a^ m and 
bi >rn of Eq.(J3J) and coefficients Q >m and r]^ m of Eq.(J2J in the following form: 

%m = Oil,rnCl,m (4) 

k,m = Pl, m T]l,m, (5) 



where scattering coefficients a and j3 are given by 

n 2 ji t m (nx) [xji >m (x)]' - ji >m (x) [nxji^m (nx)]' 
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(6) 
(7) 



Here we introduced the dimensionless parameter x defined as x = kR; f'(x) means a deriva- 
tive of / with respect to x. In what follows we will dimensionless frequency. 
The coefficients a and f3 as functions of x possess poles on a complex plane. The real parts 



of these poles give frequencies of the respective scattering resonances, while the imaginary 
parts determine the width of the resonances. It is customary to characterize the latter 
using a quality factor Q defined as Q = Kex/lmx. The resonances are characterized by 
polarization (TE or TM, depending on their origin from a or (3 coefficients, respectively), 
and two indices, /, and s, where the former is the respective angular momentum, while the 
latter is a radial quantum number enumerating poles of a and f3 with the same values of I, 
but different radial distribution of the field (the resonances are degenerate with respect to 
index m) . Following Ref 3 we will label resonances as ITEs or IT Ms. Q-factor of these 
resonances strongly depends on / reaching very high (up to 10 9 ) values for large enough 
At the same time, for modes with larger s Q decreases, sometimes by an order of magnitude. 

The resonances with high Q values correspond to so called whispering gallery modes 
(WGM), and can be sensibly described as almost stationary (long-living) normal modes of 
a sphere. The collectivization of these modes in a chain of spheres is the main focus of this 
paper. The electric field of these modes has an evanescent character in the vicinity of the 
surface of their respective spheres making the tight-binding description of the optical cou- 
pling between adjacent spheres rather accurate^. This description is based on generalizing 
the expansions for scattered wave given by Eq. (J3J) to the case of multiple spheres: 

E * = E W.mNm.K' - *i) + &U M "M( r " r *)] • ( 8 ) 

i,l,m 

where r-j is the position vector of the center of z-th sphere. The problem is, however, that 
VSH in this expansion are defined in different coordinate systems associated with the center 
of each sphere. In order to apply boundary conditions one need to rewrite them in a common 
coordinate system. This is achieved with the help of addition theore m 14 i 15 , which presents 
VSH centered at point rj in terms of VSH centered at a point rf. 

N,,„(r - r,) = Y, E K;™'fc - r<)N 1 ., m .(r - r.) + B&f (r, - r,)M,,,„,(r - r,) 

rr <»> 



l'=l m'=-V 



where A 1 ^, (rj — r^) and _B^'™,(rj — r^) are the so- called translation coefficients, which 
depends on the mutual position of the spheres. Expressions for the translation coefficients 
can be found, for instance, in RefU or Refjl2l. The latter paper gives a rather detailed 



description of methods for computing these coefficients, which is not an easy computational 
task, since for modes with large I it involves factorials of very large numbers. Using the 
addition theorem one can derive a system of equations for the expansion coefficients a\ m 
and b{ m M 



.in 



'his system of equations is valid for an arbitrary collection of spheres, 
lfll For a linear chain, in which all position vectors can be chosen 



ai,m < Cl,m + [ a i',m'K^' ~ T i) + K,m' B l',m ( T j ~ r i) \ ( 10 ) 

I j+i l',m' ) 

I j^i l',m> ) 

Indexes i and j in these equations enumerate spheres, while I, I', m, m! correspond to 
different resonances of the individual spheres. The structure of these equations can be 
understood by noting that it is equivalent to Eq. and (0), in which the incident field, 
characterized by coefficients C;,m and r]^ m is replaced by its sum with the field scattered by 
all other spheres, 
as indicated in Ref 

parallel to each other and to the axis of the chain, these equations can be simplified. If 
one chooses the Z-axis of the coordinate system along the chain, the translation coefficients 
can be shown to become J^^™ = J^[b m)m i and B\'™ = Bi5 m ^ m i. Thus, the sum over m' 
in Eq. (fTU|) and (fTTj) disappears, making the component of the angular momentum along 
the axis of the chain a conserving quantity. This fact obviously reflects the axial symmetry 
of this system. In what follows, we, for concreteness, shall assume that m = 1; results for 
other values of m can be obtained from our general formulas by recalculating parameters 
a and (3 as well as the translational coefficients. Thus, we can abridge our notations by 
dropping index m all together. At this point we also omit the terms r]^ m and Cz,m describing 
the external incident wave, which leaves us with a system of homogeneous linear equations 
and the problem of finding their normal modes instead of the problem of scattering of an 
external wave. 

Thus, we present the equations for the scattering coefficients in the following form: 

«i = «i£E [4 A (h3) + b\,B i ; (ij)] (12) 

b i = a £ E (»'. + 4>b i ; (ij)} , (is) 

where the position vectors, rj, of the spheres in the arguments of the translation coeffi- 



cients are for shortness replaced with sphere's numbers, i. Eq.()12|) and (|T^|) have been used 
by Miyazaki and Jimba in Refll2 for exact numerical analysis of a bi-sphere. They also 
developed for this system an approximate method of solution of these equations that re- 
produced results of the exact calculations with a good accuracy. Here we generalize this 
method for a system with an arbitrary, including infinite, number of spheres, where its use 
becomes crucial since exact numerical calculations grow increasingly more involved and ex- 
pensive with the increase in the number of spheres in the chain. We will show here that 
using the expanded version of the tight-binding approximation we are able to derive analyt- 
ically dispersion equations describing collective excitations of the chain, which can be solved 
numerically with minimal computational efforts. 

The reduction of Eq. (|12|) and Eq. (|13j) to an analytically tractable form is based on several 
important properties of the translation coefficients, Af, and B( . First of all, in the case of 
I 3> 1 the calculation of these coefficients can be significantly simplified with the help of 
the so-called maximum term approximation^ that allows presenting these coefficients in the 
following form 



41*,*) - -a fa + jftj _ 1} >< (l , - to* l< - i\) ("I 

flf = i^p^4d,i), (15) 

where rj defined as rj = jr^ — Tj\/R > 2 is a dimensionless distance between the centers 
of the spheres. In the case of spheres touching each other, which we shall assume in our 
numerical calculations, rj — 2. Taking into account properties of the Hankel function h[ (x) 
in Eq. flHj) ; one can see that the translation coefficients quickly decrease with the distance 
between the spheres. This property, which is a manifestation of the evanescent nature of 
the optical coupling between the spheres, allows one to keep in the sum over the spheres 
in Eq.(|12|) and (JT5j) only terms with j = i ± 1. The resulting equations constitute the 
nearest neighbors approximation for the chain of the spheres. The sum over V describes 
coupling between supermodes with different angular modes, which is the main subject of 
the present paper. These equations also contain terms proportional to translation coefficients 
B( , which are responsible for coupling between supermodes with different polarizations. It 
follows, however, from Eq. (jl5j) that for the nearest neighbors B\ -C A\ for 1,1' ^> 1. We will 
see later that the main contribution to the sum over /' comes mostly from /' < /, and the 
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cross-polarization coupling can, therefore, be neglected in most cases. In what follows we will 
also discard terms Q and rji responsible for external excitation. This leaves us with a system 
of homogeneous linear equations and the problem of finding their normal modes instead 
of the problem of scattering of an external wave. Thus, we present the final tight-binding 
equations describing TE and TM supermodes as 

1 v 

U = xxc^+ftH, (17) 

where the translation coefficients A\ = A\ (i,i + 1) play the role of the inter-site coupling 
parameters of the tight-binding approximation. They also play an important additional 
role of being coupling coefficients for supermodes with different /. These equations are very 
similar to electronic tight binding equations describing, for instance, s-p hybridizatio n 16 ' 17 . 
There is, however, a significant difference between optical and electronic problems caused 
by the non-Hermitian nature of matrix A\ describing the coupling in the optical case. We 
will see below that this peculiarity of optical tight-binding equations is responsible for the 
anomalous dispersion properties of the system under consideration. The non-Hermitian 
nature of the optical band coupling appears also in a multiple scattering (optical Korringa- 
Kohn-Rostocker) approach to this problem developed in Refill. 



III. COLLECTIVE EXCITATIONS IN THE LINEAR CHAIN OF MICRO- 
SPHERES 

A. Single band approximation 

Neglecting terms with V ^ I in equations (JTfij) and (JT7jl one obtains a simple single-band 
nearest-neighbor interaction model describing bands of the collective supermodes originating 
from ITEs or ITMs WGM's of an individual sphere with frequencies xi s . For a TE band, 
for instance, one has 

laj = A\ (a*" 1 + aj +1 ) (18) 

and the equation for TM bands has the same form with an obvious substitution of 
instead of «/. In order to isolate a single band, and convert this equation into a standard 
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FIG. 2: Positions and widths of several TE and TM supermodes obtained in a single band approx- 
imation. 

tight-binding form, we can take advantage of the fact that the expected width of the band 
of the /, s supermode, 5xi s , is rather small, 5xi s <C xi s . In this case we can expand the single 
sphere scattering parameters a or (3 around their respective resonance frequencies xi s< te or 
xis,tm- Taking into account the resonant nature of the scattering parameters, we can write 
down: 

i (x-xi s + i7i s ) ^ 



1 

Oil 



iAi 



and the similar equation for the TM polarization. Parameters A/ s in Eq. (|19J) are real valued 
positive quantities for all / and and ^i s represents the rate of the radiative decay of the 
respective WGM. Taking into account that translation parameters A\ in the maximum term 
approximation are almost purely imaginary quantities: 



-1) 



i+x 



A 



(20) 



and calculating them at the resonance frequencies xi s , one can re- write Eq.()18|) or similar 
equation for the TM polarization in the standard for the tight-binding approximation form 

(x - x ls + i lls ) a] = A ls A ls [a]" 1 + a^ 1 ) , (21) 

where Ai s = \ A\ {xi s )\. Normal modes of this system of equations are harmonic waves 

a] oc exp(iq ls Zi) (22) 

where Z{ is the z coordinate of the i-th sphere, and q\ s {x) is a Bloch wave number. It satisfies 
the dispersion equation of the same type as Eq.tJTJ), where we can now give a microscopic 



expression for the phenomenological coupling parameter k: 

k 1s = (-l) l 2A ls A ls (23) 

Dispersion equation (Eq. ([TJ) ) describes a band of excitations, which is symmetric with 
respect to the WPG frequency xi s , and whose boundaries xf are given by expression 

x^ = x ls ± k 1s (24) 

In the case of a chain consisting of a finite number of spheres, N, the dispersion equation (JT|) 
can be used to find the spectrum of the respective frequencies. To this end, one has to impose 
obvious boundary conditions for the coefficients aj, which read as a° = af +1 = 0, where we 
assume that first sphere in the chain is assigned number i — 1. These boundary conditions 
determine the allowed values of the Bloch wave number qi s : qi s r] = im/{N + 1), where n 
changes from 1 to N. Applying this result to the case of only two spheres we find that there 
are two possible values of qi s , namely qi s r) = tt/3 and qi s rj = 2n/3. The respective values 
of frequencies x]f = xi s ± Kj s /2, reproduce the results of Refill obtained for the case of a 
bi-sphere in the single mode approximation. The radiative decay characterized by 7; s makes 
the supermodes quasi-stationary; obviously the whole concept of the collective excitations 
with different q can only make sense if 7; s <C i^is- Calculations as well as experiments^ 1 ^ 
show that for modes with large I this inequality is well satisfied. 

This simple picture could describe supermodes that do not overlap spectrally since only 
in this case omitting terms with I' ^ I in Eq.(|16|) can be justified. The bands of collective 
excitations in the chain of microspheres, however, almost certainly overlap with at least one 
or more other bands. In Fig|2] we present positions and widths of a number of supermode 
bands calculated in the single band approximation. It should be understood, however, that 
because of the inherently present radiative decay, the concept of allowed and forbidden bands 
is not very well defined even for modes with relatively large Q-factors. Nevertheless, we will 
use these terms to describe spectral regions, where the wave number q would have been 
purely real or purely imaginary in the absence of the decay, respectively. In reality, the wave 
number is a complex valued quantity at all real frequencies, and a difference between allowed 
and forbidden bands manifests itself only in relative magnitudes of real and imaginary parts 
of q: Imq <C Reg within the allowed band, and Img > Reg( mod ir) within the forbidden 
bands. In the vicinity of the band boundaries, xf 1 , the real and imaginary parts are of the 
same order: Im q ~ Re q. 
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B. Inter-band coupling and band hybridization in the linear chain of spheres 



In this subsection we will take into account non-diagonal terms in Eq. (jl6j) and ()17|) . which 
are responsible for coupling between bands with different angular momentum described in 
the previous subsection. Since calculations for TE and TM modes are almost identical, we 
shall consider explicitly only TE polarization; results for TM polarization, which we will use 
for numerical examples, can be easily obtained with the help of obvious substitution of the 
parameters. 

A particular solution to the system of equations (j!6|) can be easily found in the form of 
a harmonic wave: af" = aie iqz \ where ai satisfy the following algebraic equations: 



Let us consider this equation in a frequency region occupied by a single-band supermode 
IqTEsq characterized by quantum numbers I = Iq and s = sq- We are interested in finding 
corrections to the single band dispersion law for this supermode caused by its interaction 
with other bands. In order to develop an approximate analytical solution for this problem 
we have to take into account several circumstances. First, main corrections to the diagonal 
approximation come from those I, which correspond to bands spectrally overlapping with 
IqTEsq. It should be remembered, however, that the bands discussed in the previous subsec- 
tion are broadened because of the radiative decay, and the term overlap should be understood 
with this fact in mind. One can see from FigOthat all supermodes with I > l are spectrally 
well separated from IqTEso, thus the respective terms in Eq. ()l(j)l can be considered a small 
perturbation despite the fact that A\ Q > A\ for I > l . Terms with I < l are more important 
since according to FigEI respective bands can be spectrally in close proximity to IqTEso, 
or even overlap with it. On the other hand coupling coefficients A\ Q quickly decrease when 
I becomes smaller than Z . Thus, the effect of coupling to these supermodes depends on 
an interplay between resonant enhancement due to spectral proximity, and decrease in the 
coupling parameter A\ o . 

In most practical situation, among all the bands contributing to the sum over I' there is 
just one, which we will label as l\TEsi, effecting the supermode under consideration in a 
most significant way. Usually the interaction with such a band is too strong to allow for a 
perturbative treatment. Following the terminology of Ref|l2 we will call IqTEsq and l]TEs\ 




(25) 
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significant supermodes, while the remaining ones are called bath or reservoir bands and can 
be treated perturbationally. We will show here that in the system under consideration the 
inter-band coupling can result in three qualitatively different types of the modifications of 
the dispersion laws. One, which is characteristic for weaker interaction can be called band 
renormalization, while the others, requiring a stronger coupling, are more appropriately 
called band hybridization. There might be two types of the hybridization, one arising, when 
the interacting bands are characterized by angular momentum indexes of the same parity, 
and the other one corresponding to the situation, when they have angular momentum indexes 
of different parities. We will show that the hybrid band arising in each of these cases exhibit 
qualitatively different dispersion properties. 

We will begin by developing a general theory capable of describing all possible situations, 
and then consider conditions controlling a transition from one regime to the other. As it was 
mentioned before, the interaction between IqTEsq and l{TEsi should be taken into account 
exactly, while the contribution from all other supermodes can be treated perturbatively. 
Accordingly, we present the system of Eq. (|23|) for the amplitudes of the supermodes in the 
following form: 
1 



+ 2A l ° cos {qrj) 
ai 

— + 2A l j cos {qrj) 



a i. 



+ 2 cos {qrj) A^a h + 2 cos {qrj) ^ A\ a i = 



ai 



1 + 2 cos {qrj) A^ap + 2 cos {qrj) A l h a t = 



(26) 



1 

Oil 



+ 2A\ cos {qrj) 



ai + 2 cos {qrj) A\ ay = 0. 
I'^i 



where in the first two lines we extracted from the sum over I the terms with I = Iq and 
I — lii an d wrote down the separate equations for the respective amplitudes. The last line in 
Eq. (}2l)j) represents the equation for the reservoir modes, which we solve for ai and substitute 
the solution to the first two lines of this equation. In the resulting double sum over / and /', 
we keep only terms with = l Q and /' = l±. This procedure, which constitutes the second 
order perturbation theory for the reservoir bands, results in a system of equations for the 
amplitudes of the significant supermodes, which can be presented in the following form: 



1 

— + 2A]° n cos {qrj) 
l a i 

— + 2A\\ cos {qrj) 



ai + 2 cos [qrj] A l j Q a h = 
ai x + 2 cos {qrj) A^ap = 



(27) 
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The weak interaction with the reservoir results here in the renormalization of the coupling 
coefficients according to the rule: 



A 



A = A 1 , 



V 



1 ~\~ cos {q^]\ / 



where 



A\Al 



(28) 



(29) 



is the inter-band interaction constant, and cos {q v rj) corresponds to the single band approx- 
imation for the reservoir modes defined as 



cos {q v rf) 



1 



2a v A% 



(30) 



In contrast to Sec 1111 Al we have omitted here subindex s in our notations for wave numbers 
q v . The reason for this is that there is no summation over s in any of the equations describing 
the inter-supermode coupling, such as Eq. (J23j) or (|2Hj). and the dependence on this index 
appears in an explicit form only if we expand scattering coefficients ai around a respective 
WGM. Unlike Sec JIII Al we cannot do this here, because we take into account contributions 
from bands in the frequency region, which can be relatively far away from their respective 
parent WGMs. Thus, the sum over v in Eq. (f2*Kj) contains contributions from terms for which 
the vicinities of x/ oSo and x\ xsx belong to the forbidden bands of respective supermodes, so 
that their (mostly imaginary) parameters q v , defined by Eq. ()30|) . can no longer be associated 
with a particular WGM. 

Eq. (j27|) describes new bands formed from the initial IqTEsq and l{TEsi bands. In order 
to simplify notations we will omit for now the angular momentum indices, and will label 
these new bands simply as q + and g_. We can present dispersion equations for each of q± 
in the form: 

1 



cos(q±rj) 



loh 



cos(qi r]) + cos(q h r]) ± J [cos(q k) r]) - cos(q h 7])} + 40 Wl cos(g io r/) cos(q h rj) 



(31) 



where f^ ^ is defined by Eq. (}29|) . in which we set /' = / = lo, v = li, drop the upper 
index as duplicate, and replace the coupling coefficients with their renormalized according 
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to Eq.pSjl values. Similarly, q\ are defined as in Eq.(|HUj). but also with renormalized coupling 
coefficients. 

Eq. (J31|) can be considered either as equations for wave numbers q± as functions of fre- 
quency x, or as equations for two frequencies x± as functions of wave number q. The choice 
depends on an experimental situation under consideration. In the transport experiments of 
the type carried out in Ref. BorO the frequency is fixed by an external excitation, and q is 
determined from the experiment. In this case we have to solve this equation for q± treating x 
as an independent real valued variable. The resulting wave numbers are complex, and their 
imaginary parts characterize the spatial decay of the wave along the chain. On the other 
hand, in the resonance experiments of Ref. |s| the wave number is fixed, and frequency is 
measured. In this case, we have to solve our dispersion equation for x considering q as a real 
quantity. The frequencies obtained as a result contain imaginary parts, which describe the 
spectral width of the respective resonances. It is important to understand that because the 
dispersion equation is complex, the functions Re(q±(x)) and Re(x±(q)) obtained in these two 
approaches are not inverse of each other. In this paper we will focus on finding q±(x) because 
it presents greater interest from the stand point of experiment as well as applications. 

For two special values of frequency, Eq. (j31|) can be solved exactly. Indeed, consider 
x = xi s , where cos(g/ ) is exactly equal to zero regardless of the renormalization of the 
coupling coefficients. In this case the term responsible for the inter-band coupling vanishes, 
and we obtain that one of the cosg± is also equal to zero. The same is obviously valid for 
x = xi lSl . We can conclude, thus that the centers of all bands, where q±rj = tt/2, correspond 
to frequencies of the respective parent WGM resonances regardless of the presence of the 
inter-band coupling. In resonant experiments with finite chains^ the admissible values of q 
are determined by the boundary conditions at the ends of the chain. In this situation, the 
center of the band is accessible only in systems with odd number of spheres. Therefore, 
spectra of chains with odd and even number of spheres can be distinguished from each 
other by, respectively, presence or absence of excitations at frequencies corresponding to 
single sphere whispering gallery modes. This conclusion is in complete agreement with 
observations of RefQ 

In general case we can solve Eq. (j31|) numerically by consecutive iterations. The zero itera- 
tion corresponds to neglecting the bath-induced renormalization of the coupling coefficients, 
and produces two zero-order dispersion curves q± . At the next step, q± are substituted 



14 



to Eq. (|28|) and Eq. (JHUj) producing a pair (one for q + and one for q_ ) of new values of 
the coupling coefficients and "single-band" wave numbers q\. These new values go back to 
Eq.(|31|). one to the "+" version of it, and the other to the "— " version. The procedure 
can be repeated as many times as necessary to achieve its convergence. The experience 
shows, however, that a good approximation for q± can already be obtained after only the 
first iteration. 

While the contributions from the reservoir bands can be important, the main qualitative 
characteristics of the dispersion laws of the significant supermodes can be understood from 
the zero order approximation, which is presented by Eq. (j31j) without the renormalization 
of the coupling coefficients. We are focused on the frequency region, where expression 
cos qi — cos qi x is small because this is where the main effects of the inter-band interaction 
are expected. The strongest effects occur in the vicinity of a resonance between supermodes 
lo and li, when 



if it takes place in the frequency range covering initial IqTEsq and l\TEs\ bands. It should 
be noted, however, that because of the radiative decay the band boundaries particularly 
for the band with larger s are not that very well defined, and the strong effects can take 
place even if the resonance condition is fulfilled in the forbidden band of one or even both 
of the modes, if the resonance point lies in proximity of a band boundary. It is important to 
remember also that Eq.(j321) is not equivalent to Re(g; ) = Re(^J, which is often accepted 
as a resonance condition. 

The properties of the dispersion curves in the vicinity of the resonance point are de- 
termined by two circumstances: (i) the parity of single band dispersion laws, and (ii) the 
relation between [Im [cos (q^r/)] — Im [cos (qi rj)]] 2 and cos(g; rj) cos^r/), which deter- 
mines if the band interaction is weak or strong. In most cases the overlapping supermodes 
originate from WGM with different s, and, therefore, usually, one of the imaginary parts 
in the expression above is significantly (orders of magnitude) larger than the other one. 
In this case the stronger decaying supermode determines the nature of the modification of 
the dispersion laws. It is quite obvious that the strong coupling case (band hybridization) 
would correspond to the interaction terms prevailing over the dissipative terms, and the 
weak coupling regime (band modification) would take place in the opposite situation. If 




(32) 
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the resonance condition is not fulfilled in the spectral region of interest the modification of 
the dispersion law always belongs to the weak coupling case. What is surprising, however, 
is that even in this case the modification can be so strong as to render the perturbative 
treatment of interaction between the significant modes not very accurate. 

An important peculiarity of the dispersion curves described by Eq.(|31|) distinguishing 
it from other cases of mode coupling (polaritons, for instance) is the sign of the inter- 
band interaction parameter, flj^. Using properties of the coupling coefficients A l [ it is 
easy to show that this parameter is always positive regardless of the parities of the angular 
momentum indexes I and I'. The parity of I according to Eq. (j21)) determines the sign of 
the slope of the respective dispersion laws in the vicinity of band boundaries. Taking into 
account that in the frequency region of interest cos(qi rj) and cos^q^r/) are, at least, of the 
same sign, we conclude that the entire interaction term fl^^ cos(qi rj) cos(qi 1 rj) is always 
positive regardless of the interacting bands having slopes of the same or different signs. We 
will see in the subsequent sections that this feature of the inter-band interaction, which 
results from the non-Hermitian nature of the matrix Af, is responsible for a rather unusual 
form of the dispersion laws describing hybridized bands. 

1. Weak inter-band coupling 

As an example of the band modification we consider dispersion curves emerging as a 
result of interaction between 29TM1 and 25TM2 supermodes. This choice is motivated 
by experiments of Refill, where these dispersion curves have been measured. One can see 
from FigO that the bands of these two supermodes slightly overlap. However, considering 
Re [cos (q29V)] an d Re [cos fes??)] we find that the resonance between these two bands does 
not occur, thus, we should expect in this case just a modification of the dispersion laws 
without a significant reconstruction of the band structure. 

FigEH presents the results obtained after first two iterations of numerical solutions of 
Eq.()31|) for these two modes along with the respective initial bands. The second iteration 
takes into account all reservoir modes from I = 1 to I = 40. One can see that the effect 
of the reservoir modes is almost zero in the vicinity of the center of the band, but slightly 
increases toward the edges of the initial bands, which confirms our choice of significant and 
reservoir modes. However, outside of the allowed bands of the 29TM1 and TM24, 2 modes 
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FIG. 3: Dispersion curves of 29TM1 and 24TM2 bands. Dash lines labelled 1 and 4 presents 
single-band dispersion curves, thin lines labelled 2 and 5 demonstrates dispersion cruves found 
without reservoir modes taken into account, and finally thick lines labelled 3 and 6 show dispersion 
curves with reservoir modes included 

designation of these modes as significant ones may no longer be valid because we can trespass 
to the allowed band of one or several of other modes, which we treated here as the bath. 
It is obvious that in those frequency regions the approximation used to obtain Figgis no 
longer valid. 

The shape of the modified curves can be easily understood on the basis of our general 
analysis of Eq. ljHlj) . As we explained above at qii] = tt/2 modified and initial dispersion laws 
coincide, but farther away from the center of the respective bands the modified dispersion 
laws are pushed away from the initial curves. The characteristic shape of the dispersion 
curves results from the fact that the corrections to the initial dispersion laws due to the 
inter-band coupling have different signs for 29TM1 and 25TM2 bands: its negative for the 
former and positive for the latter. The most strong modifications of both curves occur at the 
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FIG. 4: Decay rates of 29TM1 (1, 2 and 3) and 24TM2 (4, 5, and 6) bands. The labelling of 
the curves is the same as in Fig|3 The insert shows, in a magnified form, the spectral region 
corresponding to 29TM1 mode. One can see a significant increase in the radiative decay rate of 
this mode in the vicinity of the new band boundaries, which, however, decreases to the initial 
WGM value at x = X29,i 

lower frequency boundaries of their respective bands. This modification can be described 
as a shift of the boundaries of the respective bands, which is more prominent for 29TM1, 
whose initial band has much better denned boundaries because of the smaller decay. These 
boundaries are not very well defined for modified 29TM 1 band, however, which tells us about 
an increase in the radiative decay of this supermode caused by the admixture of 25TM2 
band. This increase, however, is not very dramatic, and allows keeping the concept of allowed 
and forbidden bands as a convenient tool to characterize different spectral regions. The 
frequency dependence of the decay rates, characterized by Img^ are shown in Fig0J where 
the positions of the band boundaries characterized by an abrupt change in the decay rate are 
seen much clearer. Another interesting effect seen in this figure is that the frequency 229TM1 
of the respective WGM is no longer at the center of the modified band. This is, obviously, 



18 



caused by the shift of the band boundaries, which for both low- and high-frequency edges 
occurs toward lower frequencies. This shift is also responsible for an asymmetry of the 
band, which manifests itself in the different growth of the decay rates toward left and right 
boundaries. 

We can estimate the shifts of the band boundaries for 29TM1 and 25TM2 supermodes 
if we neglect the renormalization of the coefficients A\ x , and treat the inter-band interaction 
parameter f^ as a small perturbation. Then, employing the approximation for coefficients 
oti given by Eq.([T§|). and neglecting the frequency dependence of the coupling coefficients 
we can find for new band boundaries, x^, defined as frequencies at which Re[cos(q±T])] = 1 
the following expressions 

~± _ ± -T- ^lu K ls /oo\ 

x bi — x bi + tzz r T \ ' 

1 =F cos [q v (x bi ) r)\ 



Index / in this equation refers to the band whose boundaries we are calculating, while the 
index v signifies its interacting counterpart. Since Qi u is always positive, and ni s is negative 
for both participating modes, the sign of the correction is determined by the values of cos q 
of the z/-th band at the initial boundaries of the /-th band. Since its quite likely that these 
boundaries lie in the forbidden band of the v-th band, the values of these cosines are not 
limited by unity. In the case of 29TM1 supermode interacting with 25TM2 we find that 

cos [225(^29,1)77] ~ 0- 18 cos [fe(^29,i)^] ~ -1-46 

cos [q 2 9(xt 5:2 )v] » 2.53 cos [fefe^) 7 ?] » -0.5174 

These calculations show that the corrections to the upper and lower boundaries are both 
negative for 29TM1 band, and both positive for 25TM2 band in complete agreement with 
numerical calculations. 

As it was mentioned we choose the supermodes 29TM1 and 25TM2 for illustration of our 
general results mostly because their dispersion curves were observed experimentally in Reffl 
However, before comparing the experimental results with our theory we should recall that 
the experiment in Ref Q was conducted in a resonance configuration, when qi was fixed by the 
conditions of the experiment, and real and imaginary parts of frequency were measured via 
positions and widths of the respective resonances. This experimental setup is different from 
what was assumed in our theoretical analysis. Therefore, the comparison makes sense only 
away from the band boundaries, where decay rates are relatively small and the inverse of 
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FIG. 5: Experimenta 
of the figures in Ref|£ 



dispersion curves of 29TM1 and TM24, 2 bands obtained by digitizing one 



the experimentally observed quantity xi(q) is close to qi(x) studied here. Fig demonstrates 
the curves q29,i(x) an d 925,2 {%) produced by digitizing Fig.2b from RefQ. Comparing this 
figure with results of our calculations, FiglBlone can see an excellent qualitative agreement 
between the theory and the experiment. 



2. Strong inter-band coupling: band hybridization 

In order to illustrate the effects of the strong inter-band coupling we consider the inter- 
action between modes 39TE1 and 3ATE2. The resonance condition (|32p for these modes is 
fulfilled at frequency x r , which is within allowed spectral region for both bands, albeit rather 
close to their respective band boundaries. Therefore, we should expect the manifestation 
of the strong coupling effects in this case. Another important feature distinguishing this 
pair of modes is the different parity of their angular momentum indexes. As a result, the 
initial (single-band) dispersion curves of these supermodes are characterized in the vicinity 
of the resonance point by slopes of opposite signs. In the case of regular interacting modes 
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such as phonon- or exciton-polaritons, one would expect in a situation like this a normal 
anti-crossing behavior resulting in the opening of the polariton band-gap in the spectrum of 
the system. This would have happened in the case considered here as well, if the interaction 
term in Eq. (j31|) were negative. It, however, is always positive as it was explained above 
and, therefore, no spectral gap arises. The shape of the dispersion curves emerging in this 
situation and shown in Fig|B] is rather unusual. The parts of the initial dispersion curves 
located below the crossing point are pushed downward from it and connect in a continuous 
manner to form a new hybrid band. The remaining parts of the initial dispersion laws are 
pushed upward such that the wave numbers corresponding to the resultant second dispersion 
curve become complex valued in the entire spectral region under consideration, with its real 
part remaining equal to n for all considered frequencies. In other words, the spectral region 
covering initial 39TE1 and 3ATE2 bands is a forbidden band for the second dispersion curve. 
Thus, the inter-band coupling turns two co-existing initial bands into a single hybrid band 
with rather unusual dispersion properties. Instead of a spectral gap we have here a region 
of forbidden values of wave numbers, which extends from the maximum value of the wave 
number on the lower dispersion curve, q max , all the way to the boundary of the Brillouin 
zone. This means that the states characterized by the wave numbers from the forbidden 
region cannot have a real valued frequency, i.e. do not correspond to stationary solutions of 
the respective dynamic equations. 

The new allowed hybrid band is a combination of the initial bands for the two interacting 
supermodes and covers a spectral region previously occupied by 39T_E1 and 34TE2. The 
shifts of the initial band boundaries (the low frequency boundary for 34TE2 and higher fre- 
quency boundary for 39TE1) can be estimated with the help of the same expression (}3"3"j) that 
we obtained for the weak coupling case. This is possible because these band boundaries lie 
far away from the resonance point and we can use perturbation theory. Taking into account 
that the sign of parameters k is now different for the two supermodes we find that the the 
low frequency boundary of 3ATE2 shifts toward higher frequencies by approximately 0.025, 
while the higher frequency boundary of 39TE1 supermode shifts toward lower frequencies 
by approximately 0.01. These numbers are in good agreement with numerical calculations 
shown in FiglHl This figure also shows that taking into account the reservoir modes (second 
iteration of our procedure) does not change properties of the hybrid mode too much. 

While the existence of the wave number gap is the general property of the dispersion 
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FIG. 6: Dispersion laws describing hybridized bands arising as a result of interaction between 
39TE1 and 2>ATE2 supermodes. Curves 1 and 2 represents initial unperturbed dispersion laws, 3 
and 5 are dispersion curves obtained without accounting for interaction with reservoir modes, and 
curves 4 and 6 depict dispersion laws of two branches with reservoir modes taken into account. 
The upper branch is almost completely pushed to the forbidden region. 

equation (|31|) for the significant supermodes characterized by angular momentum index of 
different parity, the complete vanishing of the upper hybrid band is specific for the pairs of 
modes under consideration here. In order to understand the properties of these bands better 
we can estimate the distance between the lower and upper branches exactly at the resonance 
point neglecting the effects due to the reservoir modes. This distance, E, determines the 
position of the maximum q ma x for the lower branch, and the minimum q m in for the upper 
branch. It is more convenient to work with cos q± than with the wave numbers themselves, 
so we will define E via equation: 

cos (q m axV) = c os (q r rj) + E 

(34) 

cos (q m inV) = cos (q r rj) - E 
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The magnitude of X can be estimated from Eq. (|31|) , which gives 



Re 



^4Q hlo - r 2 - 4^ /l/o r 



where we took into account that the resonance occurs close to the band boundaries, so that 
Re[cosg; ] ~ Re [cos 9^ ] ~ —1, and introduced T = Imfcos^J, which represents the super- 
mode with the largest radiative decay (34TE2 in the case under consideration). However, 
the radiative decay of the 34TE2 supermode is still so small that T 2 <C in this case, 
and £ can be approximated as 

S = (36) 

In the particular example of 39TE1 and 3ATE2 bands the resonance occurs at a point, where 
cosg r is so close to the boundary value of —1, that cos (q m inV) becomes less than —1 making 
the respective wave numbers for the entire frequency region imaginary. Should the crossing 
point of the two interacting supermodes lie farther away from the band boundaries, the part 
of the upper hybrid band could also survive, but it would be restricted by the frequency 
region in the vicinity of the resonance frequency. The cos (q m axV) on the other hand is 
pushed farther away from the boundary, so that the entire lower band remains allowed, 
and the respective maximum allowed value of q can be estimated as q ma xV ~ qrV — \/2£. 
Comparison of these estimates with numerical results shows that they give a relatively good 
approximation for the respective quantities in the case under consideration. 

One of the consequences of having a gap for wave numbers, is that formally speaking, the 
group velocity of the excitation described by the hybrid dispersion law diverges at q = q ma x- 
This form of a dispersion curve, therefore, raises a question if it is consistent with casuality. 
This question was discussed in Refllfli where a similar form of the dispersion curve was 
found in an one-dimensional resonant photonic crystal. It was shown in that work that the 
wave-number gap and accompanying it infinite group velocity do not contradict causality, 
if attenuation of the respective excitation is properly accounted for. Similar anomalous 
behavior of the group velocity (without the wave number gap, of course) has been known 
for a long time, for instance, in the case of exciton-polaritons. 2 - A more general statement 
that abnormal behavior of group velocity should be expected in the regions of resonant 
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absorption in all dispersive dielectrics was proven in Ref 

In spite of formal analogy between the anomalous dispersion law presented in FiglHl and 
the results of Refill these two situations significantly differ from each other. The wave 
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number gap found in Refll9l is caused by the frequency dispersion of the respective dielectric 
medium, and coincides with regions of anomalous dispersion and resonant absorption. The 
origin of the gap in our case is completely different and can be traced to the combination 
of two factors: different signs of the slopes of the initial dispersion curves and positive 
value of the respective coupling constant. The latter is caused, as we already discussed, 
by a non-Hermitian nature of the matrix describing the inter-band coupling in the case 
under consideration. Accordingly, the decay rate of the photonic supermodes considered 
here does not show any resonant enhancement in the vicinity of the wave-number gap. We 
present the frequency dependence of this quantity, which remains surprisingly small in all 
frequency range corresponding to the allowed hybrid band, and demonstrates a weak non- 
monotonic dependence on frequency with a flat minimum (FigEJ). We can explain this 
behavior qualitatively by reminding that at the centers of both initial bands the interaction 
between them vanishes, so the decay rate at these two points coincides with the decay rates 
of the original WGM resonances. The crossing of the curves 2,3 and 4 in the insert in Fig0 
corresponds to the initial 34TE2 WGM, and the second (counting from the left) minimum 
on curve 4 corresponds to x^qtei- The first minimum on this curve appears only after the 
interaction with reservoir modes is taken into account and presents another peculiarity of 
the system under consideration. Normally, one would expect that the reservoir bands can 
only increase the radiative rate, while here we observe an opposite behavior. A possible 
qualitative explanation of this effect can be offered on the basis of the following arguments. 
Let us assume, for an instance, that instead of considering the experimental situation, in 
which frequency is considered real and the wave number complex, we deal with the resonance 
type of experiment with real q and complex x. In this case it is reasonable to assume, and 
the results of Refill support this assumption, that the the reservoir modes would play its 
regular role and increase the width of the respective resonances, Ax. The transition from this 
description to the one used in this paper, with complex valued q can then be approximately 
carried out with the help of the following expression: Img ~ (dRe q / dx) Ax . Thus, the 
frequency regions with flat dispersion curve should be characterized by decreasing imaginary 
part of the wave number. In order to understand why this minimum appears only after 
reservoir bands are taken into account, one needs to compare respective dispersion curves. 
From FigEJ one can see that without the reservoir modes the maximum of the dispersion 
curve is relatively narrow. The reservoir modes make the maximum wider, which means 
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FIG. 7: Radiative decay rates for the allowed hybridized band. The labeling of the lines is the 
same as in FigEl The insert magnifies the region around the center of 39TE1 band. 

that a broader frequency interval is characterized by small values of dKeq/dx. This causes 
a faster decrease of the spatial decay rate when frequency shifts toward the allowed region 
for former 39TE1 band, which can be seen in Fig[7| This tendency reverses for frequencies 
greater than x r , when the group velocity starts decreasing. When, however, the frequency 
approaches x 3 gxEi, the decay rate decreases again resulting in a curve with two minima, one 
at x r , and the other one at x^tei- 



IV. CONCLUSION 



In this paper we have developed a theory of the inter-band coupling in the system of 
dielectric spheres forming a one-dimensional chain. Our objective was to obtain analytical 
expressions allowing one to study band structure and dispersion properties of the collective 
excitations of this system with inter-band mixing effects taken into account. We developed 
an approximation scheme generalizing an approach of Ref jjljji which was originally applied to 
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the case of two coupled spheres and derived a general dispersion equation describing the band 
structure emerging as a result of the inter-band interaction. We showed that there might 
exist three qualitatively distinct regimes of such coupling depending upon the properties of 
the initial single-band dispersion curves of the interacting supermodes, and the strength of 
the interaction. Our general results were applied to two particular examples illustrating two 
of the possible manifestations of the inter-mode coupling. 

As an illustration of the regime of weak band modification we considered frequency region 
corresponding to initial supermodes TE29, 1 and 25TM2, whose initial dispersion curves do 
not cross each other. Nevertheless, we showed that the inter-band coupling significantly 
modifies the dispersion curves of these supermodes and is responsible for their characteristic 
shape observed in experiments of Ref|a What is interesting in this example is that the 
modification of the dispersion laws mimics the anti-crossing behavior typical for resonantly 
interacting excitations. At the same time, as it was mentioned above no crossing resonance 
takes place in this case. We explained the shape of these curves as a result of two effects. 
First, we showed that the position of frequency corresponding to the center of the band, 
where wave number q = ir/(2rj), where 77 is dimensionless period of the structure, is pinned 
to the frequency of the parent WGM mode corresponding to a given band. Second, we 
demonstrated that the inter-band coupling shifts the boundary of the respective bands in 
the opposite directions. As a result, one of the dispersion curves bends upward, while 
the other one bends downward imitating the anti-crossing behavior. The results of the 
numerical calculations of these dispersion curves show very good qualitative agreement with 
experimental results of Reffl 

More interesting situation was found, however, in the spectral region corresponding to 
supermodes 39TE1 and 34TE2. These modes exhibit true crossing resonance very close to 
the boundaries of the both bands. It is also important that in the vicinity of the crossing 
point their initial dispersion curves have slopes of opposite signs. The interaction between 
these supermodes gives rise to a new hybrid band with a highly unusual dispersion charac- 
teristics. Instead of a standard avoiding crossing kind of behavior accompanied by opening 
of a gap in the frequency spectrum of the system, we found a gap in the allowed values of 
the wave numbers. More specifically, instead of two initial crossing dispersion curves, there 
emerges one curve characterized by a non-monotonic dependence of the wave number of 
frequency. The maximum value of the wave number is significantly lower than vr/r/, which 
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determines traditional band boundaries. Thus, there is a range of wave numbers that do not 
correspond to excitations with any real valued frequency. Such a band structure is not re- 
lated to the presence of the radiative decay, which remains very small throughout the entire 
band of frequencies corresponding to this hybrid band. We traced the origin of this effect to 
the non-Hermitian nature of the inter-band coupling typical for electromagnetic problems. 
The phenomenon of the wave number band-gap could be observed experimentally in an 
experiment similar to one carried out in Refj8l In that experiment the fluorescent spectra 
were observed for chains with varying number of spheres. In normal situation, adding a 
sphere results in appearance of an additional peak on the spectrum so that each peak can 
be identified with a respective wave number. The wave number gap would manifest itself 
in this experiment as a failure of a new peak to appear after addition of a sphere to the 
chain. The dependence of the fluorescence spectrum upon the number of spheres in this 
situation is a non-trivial problem requiring a separate consideration. Other possible exper- 
imental manifestations of the predicted band structure may include unusual behavior of a 
pulse propagating along the chain, and angular dependence of the radiation emitted by the 
chain. More detailed analysis of these effects is also outside of the scope of this paper and 
will be carried out in the future. 
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